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ABSTRACT 

How rapidly collapsing parsec-scale massive molecular clumps feed high-mass stars, and how they 
fragment to form OB clusters, have been outstanding questions in the field of star-formation. In this 
work, we report the resolved structures and kinematics of the approximately face-on, rotating massive 
molecular clump, G33.92-f-0.11. Our high resolution Atacama Large Millimeter/submillimeter Array 
(ALMA) images show that the spiral arm-like gas overdensities form in the eccentric gas accretion 
streams. Eirst, we resolved that the dominant part of the ^0.6 pc scale massive molecular clump 
(3.01^‘4'10^ Mq) G33.92+0.11 a is tangled with several 0.5-1 pc size molecular arms spiraling around 
it, which may be connected further to exterior gas accretion streams. Within G33.92-f0.llA, we 
resolved the ^0.1 pc width gas mini-arms connecting with the two central massive (100-300 Mq) 
molecular cores. The kinematics of arms and cores elucidate a coherent accretion flow continuing 
from large to small scales. We demonstrate that the large molecular arms are indeed the cradles of 
dense cores, which are likely current or future sites of high-mass star formation. Since these deeply 
embedded massive molecular clumps preferentially form the highest mass stars in the clusters, we 
argue that dense cores fed by or formed within molecular arms play a key role in making the upper 
end of the stellar and core mass functions. 

Subject headings: evolution-ISM: individual objects (G33.92+0.11)-stars: formation 


1. INTRODUGTION 

Magnetic field, turbulence, and (proto)stellar feedback 
are well-known physical mechanisms that regulate molec¬ 
ular cloud fragmentation, and play a crucial role in de¬ 
termining star-forming efficiency and steiiar/core mass 
function (Wang et ai. 2010; Ghen & Ostriker 2014; Tan 
et ai. 2014; Van Loo et ai. 2014). Recentiy, much atten¬ 
tion is drawn to the hierarchicai cioud contraction regu- 
iated by fiiamentary or sheet-iike gas structures (Burkert 
& Hartmann 2004; Zhang et ai. 2009; Myers 2009, 2011; 
Schneider et ai. 2010; Liu et ai. 2012a; Toaia et ai. 2012; 
Busquet et ai. 2013; Dobbs et ai. 2013). 

On the other hand, for (giant) moiecuiar ciouds that 
are undergoing rapid giobai coliapse (e.g. Liu et ai. 
2012a; Hartmann et ai. 2012; Vazquez-Semadeni et ai. 
2007; Howard et ai. 2014), the amphfied effects of the 
initiai specific anguiar momentum may iead to a distinc¬ 
tive mode of fragmentation in cioud centers (Keto et ai. 
1991). The rotationaiiy supported structures can con¬ 
tinue accumuiating gas, and can have sufficient time to 
fragment. The arm-hke gas iocai over densities, which is 
conducive to fragmentation, may aiso be shock induced 
at where the eccentric accretion streams coiiide with each 
other. Indeed, the observed marginaiiy centrifugaiiy sup¬ 
ported accretion flows (e.g. Gaivan-Madrid et ai. 2009; 
Liu et ai. 2010; Liu et ai. 2013; Gesaroni et ai. 2011) 
have argued that the gas dynamics can be dictated by an¬ 
guiar momentum in the centrai <1 pc regions of moiec¬ 


uiar ciouds. Sever ai interesting cases have shown frag¬ 
mentation in rapidiy rotation systems. Eor exampie, the 
spatiaiiy resoived ciumpy gas toroid in the center of the 
iow-mass ciuster-forming region LI287 and the center of 
the intermediate-mass star-forming region NGG6334V 
(Juarez et ai. in prep.), and the ring-iike distribution of 
the uitra compact (UG) Hll regions (Ghurchweii 2002) 
in the centrai ^2 pc of the Gaiactic mini-starburst re¬ 
gion W49N (Weich et ai. 1987; De Free et ai. 1997; 
Gaivan-Madrid et ai. 2013). Erom the observationai 
point of view, two outstanding questions remain: a) How 
the iarge scaie streams converge onto these rotating sys¬ 
tems, and b) how the dense cores form subsequentiy. 

Since the spinning massive moiecuiar ciumps are geo- 
metricaiiy flattened, the morphoiogicaiiy weii-resoivabie 
cases are iimited to those which are far from an edge-on 
projection. Their detaiied gas kinematics is oniy now ac- 
cessibie by high spectrai resoiution observations of moiec¬ 
uiar iine, thanks to the unprecedented sensitivity of the 
Atacama Large Miiiimeter Array (ALMA). We seiected 
to observe the I/5oZ^2.5xlO^ Lq OB ciuster forming re¬ 
gion G33.92-^0.11, which is at a distance of kpc 

(Eish et ai. 2003)^. Despite the strong free-free emission 
associated with the aiready existing embedded OB cius¬ 
ter, our previous Submiiiimeter Array (SMA) chemicai 
studies suggested that this source is at an eariy evoiu- 

^ The visr of this source is close to the tangent point, and there¬ 
fore does not have near-far ambiguity in the kinematics distance. 
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Fig. 1. — The SHARC2/CSO+SPIRE/Herschel 0.35 mm continuum image (^hPBW= 9".6). The significant contours are not presented 
with identical intervals, for the sake of tracing the dust/gas column density. Black contours are 250 mJybeam”^ (5cr)x[4.5, 6.0, 7.5, 9.0, 
12.0, 15.0, 18.0, 21.0, 27.0, 33.0]. White contours are 250 mJybeam”^ (5f7)x[40, 60, 80, 100, 150, 200]. The circle shows the ALMA 12m 
array primary beamwidth at 15% power (~37" in diameter). We refer to the two regions enclosed by the green dashed lines as G33-north 
and G33-south. Golor bar is in units of Jybearn”^. 


tionary stage, and may not yet be seriously disturbed 
by (proto)stellar feedback (Liu et al. 2012b). In addi¬ 
tion, previous lower angular resolution Berkeley Illinois 
Maryland Array (BIMA), Very Large Array (VLA), and 
SMA observations of the molecular line emission consis¬ 
tently found that the measured molecular gas mass is 10 
times larger than the virial mass (Watt & Mundy 1999; 
Liu et al. 2012b). This can be naturally interpreted by 
the nearly face-on projection. Therefore we expect min¬ 
imized confusion among structures in the line-of-sight. 

The major improvement of the present work over pre¬ 
vious observations (e.g. Liu et al. 2012b), is the ^8 times 
better spectral resolution ( r\j 0.17 kms ^) in the CH 3 CN 
J=12-ll observations. In addition, the unprecedented 
sensitivity of ALMA allows imaging various optically 
thiner dense molecular gas tracer at sub-arcsecond reso¬ 
lution. This permits discriminating the blended molec¬ 


ular gas arms^, which is important for diagnosing the 
details of molecular gas kinematics. In this paper, we 
focus on the 1.3 mm dust continuum emission, which is 
a reliable tracer of the molecular gas column density; the 
DCN 3-2 line, which was observed to be less abundant 
in hot molecular cores (Rodgers & Millar 1996; Hatchell 
et al. 1998); the 5-4 line which is only excited in 
warm and dense gas and may be shock-enhanced; and the 

^ We follow the existing nomenclature in the literature (e.g., 
Zhang et al. 2009; Liu et al. 2012ab). In this way. Massive 
molecular clumps refer to structures with sizes of ~0.5-l pc, mas¬ 
sive molecular cores refer to the <0.1 pc size structures embedded 
within a clump, and condensations refer to the distinct molecular 
substructures within a core. Fragmentation refers to the dynami¬ 
cal process that produces or enhances multiplicity. Molecular fil¬ 
aments refer to the geometrically elongated molecular structures, 
and molecular arms refer to segments of molecular filaments that 
are located within the <1 pc radii of molecular clumps and may 
not be fully embedded within molecular clumps. 
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Fig. 2.— The free-free model-subtracted dust coutiuuum image, takeu by the ALMA 12m+ACA array observatious 
(^maj P.A.=90°). Coutour levels are 0.6 mJybeam“^ (3a) x[-l, 1, 4, 8, 12, 24, 48, 60]. The peak emissiou iu 

this image is 46.6 mJybeam“^ at core A2. The color image does uot preseut the full iuteusity rauge because the very bright emissiou 
regious are compact. We label the source C iu the east, aud label the (A) cores discussed iu the text with uumbers. Arrows iudicate the 
molecular arms. The color bar is iu uuits of mJybeam”^. 


hot molecular core tracer, the CH 3 CN J= 12 -ll K-ladders 
(Sutton et al. 1986; MacKay 1999; Araya et al. 2005). 
Based on these observations, we are able to further differ¬ 
entiate the kinematics gas streams at different physical 
conditions, which may be well blended in the spatial or 
velocity domain. Our ALMA observations also cover the 
recombination lines H30(a and He30(a, which trace the 
ionized gas. We additionally present the SHARC 2 /CSO 
(Dowell et al. 2003) and SPIRE/Herschel^ (Griffin et al. 
2010) observations of the 350 jam continuum emission, to 
outline the overall molecular cloud geometry. Details of 
our observations are provided in Section [2l We present 
the obtained continuum image and the molecular line 
data cubes in Section [3l Section I4d1 and will address 
the observed matter distribution and the gravitational 

^ Herschel is au ESA space observatory with scieuce iustrumeuts 
provided by Europeau-led Priucipal luvestigator cousortia aud with 
importaut participatiou from NASA. 


stability of the resolved gas structures. Our interpreta¬ 
tion of the observed velocity field is provided in Section 
14.31 The phys ical implication of our results is discussed 
in Section [T4l A brief summary will be given in Section 
[5l The identified line species from our ALMA observa¬ 
tions are provided in Appendix A. 

Based on our IRAM-30m mapping observations of the 
^^GO 2-1 and G^^O 2-1 line (Liu et al. 2012c), we are 
sure that all gas structures traced by the dust continuum 
emission addressed in this manuscript are physically as¬ 
sociated. We will present the larger scale gas kinematics 
traced by the GO 2-1 isotopologues, the details of core 
mass distribution, the ionized gas kinematics traced by 
the H30(a and He30a lines, as well as the chemistry, in 
separated papers. 
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Fig. 3.— The spectra taken with the ALMA 12m+ACA Array (Section l2.ID and the identified spectral lines (see also Table [T|. Due to 
bandpass defects, the 12m Array data were lost in 300 spectral channels in the 231 GHz spectral window. 


2. OBSERVATIONS 

2 . 1 . ALMA 12m+ACA Array 1.3 mm 

The ALMA 12 m Array (i.e. 12 m dish size) obser¬ 

vations were carried out on 2014 May 04, with ^36 
good antennas. The pointing and phase referencing cen¬ 
ter is R.A. (J2000) =18^52“50®.272, and decl. (J2000) 
=00°55'29".604. We used two 234.4 MHz wide spectral 
windows (channel spacing 61 kHz, ^0.085 kms“^), and 
two 1875.0 MHz wide spectral windows (channel spacing 
488 kHz, ^0.65 kms“^), which tracked the velocity of 
i;/sr^l07.6 kms“^ on our target source. 

The central frequency of these spectral windows 
are 231.220690 GHz (^^CS 5-4), 231.900928 GHz 
(H30(a), 220.679320 GHz (GH 3 CN J= 12 -ll, K=4), and 
217.104980 GHz (SiO 5-4), respectively. The uv sam¬ 
pling range of these observations is 13-430 kA. The over¬ 
all on-source time was 43.7 minutes, and the T^^^ was 
in range of 60-150 K. We observed Titan, J1851+0035, 
and J1751+0939, for absolute flux, gain, and passband 
calibrations, respectively. 

The ALMA Gompact Array (AGA; 7m dishes) observa¬ 
tions were carried out on 2014 May 03 and 2014 May 04, 
with ^10 available antennas^. The pointing and phase 
referencing center, and the spectral setup, are the same 
with the 12 m array observations. These AGA observa- 

^ Not all antennas were always valid for the observations because 
of the antenna shadowing. 


tions cover a uv sampling range of 7.2-47 kA. The overall 
on-source time is 24.4 minutes on May 03, and 73.2 min¬ 
utes on May 04. The T^^^ was in range of 60-150 K. 
We observed Titan/Neptune, J1851+0035, and J1733- 
1304/J1924-2914, for absolute flux, gain, and passband 
calibrations, respectively. 

The data were calibrated, phase self-calibrated (2 itera¬ 
tions), and jointly imaged using the GAS A package (Mc- 
Mullin et al. 2007) version 4.2.0 (release r28322). In the 
end, we did not implement the phase self-calibration solu¬ 
tion to the 12 m Array data because there was no obvious 
improvement in image quality. We did not self-calibrate 
the visibility amplitude because of a significant missing 
flux (>24%) as compared with previous SMA observa¬ 
tions in the same region. The continuum image was con¬ 
structed by jointly imaging all spectral line-free channels 
using the multi-frequency imaging (mf s) algorithm. The 
Briggs Robust=0 weighted 12 m+ACA continuum im¬ 
age has a synthesized beam of ^maj ><^min=0''.67x0".47 
(P.A.=81°), and the dynamic range limited rms noise 
level of 0.2 mJybeam”^. We have corrected the primary 
beam attenuation, and therefore the images are noisier 
in outer region. 

The selected ^^GS 5-4 line was simultaneously cov¬ 
ered in one 234.4 MHz spectral window, and a 1875.0 
MHz spectral window. Due to a correlator failure, we 
lost the ^^GS 5-4 data in the 234.4 MHz spectral win¬ 
dow. The GH 3 GN J= 12 -ll K-ladders were observed in 
a 234.4 MHz spectral window; the DGN 3-2 and SiO 5- 
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Fig. 4.— The velocity integrated intensity maps of the selected molecular lines (color), overlaid on the 1.3 mm dust continuum image 
(contours). Contours are 0.6 mjybeam”^ (3cr)x[l, 4, 12, 24, 48, 60] in the and DCN line panels, and are 0.6 mjy beam“^ x [4, 12, 

24, 48, 60] in the panels for the CH3CN K-ladders. 


4 lines were observed in 1875.0 MHz spectral windows. 
We performed Briggs Robust =2 weighted 12 m+ACA 
spectral line imaging to enhance the sensitivity to dif¬ 
fuse emission, which yielded the synthesized beam of 
^maj x^min=0".78x0".67. The standard 1-iteration spec¬ 
tral banning smoothing procedure, which can better sup¬ 
press the passband ripple, will yield ^0.17 kms“^ ve¬ 
locity resolution, and ^1.3 kms“^ velocity resolutions 
for the 234.4 MHz and the 1875.0 MHz spectral win¬ 
dows, respectively. However, we imaged the 5-4 

and the DCN 3-2 lines with ^0.68 kms“^ channel spac¬ 
ing, intending to better present the kinematics infor¬ 
mation. The CH 3 CN J= 12 -ll K-ladders were imaged 
with the 0.17 kms“^ velocity resolution. We achieved 
the rms noise level of 1.6 mJybeam”^ (70 mK) and 3.5 
mJybeam”^ (170 mK) in the ^0.65 kms“^ and 0.17 
kms“^ velocity channels, respectively. There is no zero¬ 
spacing information available for our 1.3 mm continuum 
map. We have compared our ALMA observations with 
the existing SMA observations of the 1.3 mm continuum 
emission, the CSO SHARC 2 bolometric single-dish ob¬ 
servations of the 0.35 mm continuum emission (see be¬ 
low), and the VLA observations of the NH 3 lines (Liu et 
al. 2012b). We confirmed that the resolved morphology 


by ALMA is not significantly biased by the missing flux 
on the ^30" scale. The ALMA 12 m+ACA Array obser¬ 
vations of the dense molecular gas tracers are less biased 
by missing zero spacing because the emission regions are 
less extended in the individual velocity channels. Pre¬ 
sentation of the ALMA images will be limited to the 
15% power of the primary beam width, if not specifically 
mentioned. 


2 . 2 . Herschel SPIRE 0.35 mm 

For this study we made use of a 350 /im map of 
G33.92-b0.11 obtained by the Herschel space telescope. 
We used the level 3 (science grade) processed mosaics of 
the SPIRE-PAGS paralell observations of a much larger 
map (observation ID 1342207031) of the W49 complex 
observed in October 2010. No significant glitches are 
visible in the region of interest. Level 3 maps are pro¬ 
duced from a combination of all available contiguous ob¬ 
servations of a program. The Level 3 mosaics have opti¬ 
mal WGS solutions and analysis quality grade with zero 
points calibrations estimated using Planck Obsevatory 
maps, and archived directly in units of surface bright¬ 
ness (MJy/sr). We applied, however, a simple astromet- 










































6 


Liu et al. 


ric check by re-creating the mosaic with Montage. and 
comparing positions of compact features against a WISE 
Observatory 3.6 micron map. 


2.3. CSO SHARC2 0.35 mm 

High angular resolution continuum observations at 
0.35 mm were carried out using the SHARC 2 bolometer 
array, installed on the Caltech Submillimeter Observa¬ 
tory (CSO). The array consists of 12x32 pixels (approx¬ 
imately 85% of these pixels work well). The simultaneous 
field of view (FOV) provided by this array is 2'.59x0'.97, 
and the diffraction limited beam size is ^ 8 ". 8 . 

The data were acquired on 2014 March 24 
(^^225 ghz ^0.06). The telescope pointing and focus¬ 
ing were checked every 1.5-2.5 hours. Mars was observed 
for absolute flux calibration. We used the standard 
10 ' xlO' on-the-fly (OTF) box scanning pattern, centered 
on R.A. (J2000) =18^52^50^.272, and deck (J2000) 
=00°55'29".604. The total on-source time was 30 
minutes. Data calibration was performed using the 
CRUSH software package (Kovacs 2008). In addition, 
we used the MIRIAD (Sault et al. 1995) task immerge 
to perform the weighted sum of the SHARC 2 0.35 mm 
image with the Herschel SPIRE image (see Section 
1221) at the similar wavelength, in the fourier domain. 
This minimizes the missing flux in the ground based 
single-dish bolometric observations due to the sky 
subtraction, which typically causes defects on angular 
scales larger than the simultaneous field of view (or the 
angular throw in the nodding observations). The final 
map was smoothed to an angular resolution of 9". 6 , and 
the rms noise level achieved was ^60 mJybeam”^. In 
addition, we aligned the final 0.35 mm image with the 
Herschel, ALMA, and SMA images. 


3. RESULTS 

3.1. Continuum emission 

The gas mass is usually estimated with dust thermal 
continuum emission. While the continuum emission at 
the 0.35 mm wavelength is dominated by dust, the con¬ 
tinuum emission at 1.3 mm can have contributions from 
both the dust thermal emission and the free-free emission 
from the ionized gas. To separate these two components, 
we first estimate the free-free emission based on the ob¬ 
servations of the hydrogen recombination line H30 g line 
(see Section [321) • Assuming that the free-free continuum 
emission at 1.3 mm is optically thin, the peak hydrogen 
radio recombination line to free-free continuum flux ratio 
(Tl/Tc) is given by 




, where 5v is the FWHM of the hydrogen radio recom¬ 
bination line, n is the observing frequency, and Tq is the 
electron temperature. We adopt the nominal value of 
Te^8000 K for the electron temperature, and estimate 
Tl and 5v based on the observed peak brightness and 
the second moment in the H30 g line image cube. The 
assumed Tg is an upper limit, since the continuum is over 
subtracted on small scale when using higher values for Tg. 


^ http://montage.ipac.caltech.edu 



Fig. 5. — Top: The color composited image generated from the 
integrated intensity maps of DON 3-2 (magenta), 5-4 (cyan), 

and CH 3 CN J=12-ll K=0 line (yellow). Bottom: The dust con¬ 
tinuum image taken by the ALMA 12 m+ACA array observations. 
These two panels are presented on the same spatial scale. 


The resultant 1.3 mm free-free continuum emission model 
has a peak flux density of 16.7 mJybeam”^. This free- 
free emission model is then subtracted from the 1.3 mm 
continuum image to yield the 1.3 mm dust continuum 
emission image. We smoothed the 1.3 mm continuum 
image to the angular resolution of the H30 g line image 
before subtracting the free-free emission model. 

Figured] and [2] shows the high resolution 0.35 mm con¬ 
tinuum image and the 1.3 mm dust continuum image. 
Both continuum images trace the overall gas column 
density, but the later is optically thinner (i.e. bright- 
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Fig. 6. — The intensity weighted average velocity maps of the selected molecular lines (color), overlaid on the 1.3 mm dust continuum 
image (contours). Contours are the same as those in Figured The velocities are relative to 77^5^=107.6 kms“^. We note that the blueshifted 
CH 3 CN K=0 component may be blended with the redshifted K=1 component (see Figure O. 


ness temperature Tg^T^ust) thus can penetrate inside 
the densest regions. The 0.35 mm continuum image 
peaks at the location of the ^10 ^ Mq massive molec¬ 
ular gas clump G33.92+0.11 A (Watt & Mundy 1999; 
Liu et al. 2012b). The peak 0.35 mm flux density is ^73 
Jybearn”^ (^1.3 K, but can be beam diluted). The 0.35 
mm image resolves the massive companion G33.92-h0.11 
B northeast of G33.92+0.11 A. The dominant large scale 
features are the projected ^5 pc scale elongated struc¬ 
tures, connecting from northeast and southwest (G33- 
north and G33-south, hereafter) to G33.92+0.11 A and 
B. There might be gas filaments connecting G33.92+0.1 
A from the west, however, the individual components 
are not as clearly separated because of confusion with 
foreground/background gas structures. 

The 0.6" resolution 1.3 mm dust continuum image 
(Figure [2]) resolved abundant structures within the in¬ 
ner parsec scale radius. The most prominent ones 
are the two 100-300 Mq massive cores Al and A2 
located at the center, and the five molecular arms 
roughly following the paths of Al—>A2—>A3 (arm-cl), 
mini-arm-A2w—>A4^A5 (arm-c2; see Figure 01 [6]), 
A4^A6^A11^A12->A13 (arm-c3), A1^A9 (arm-c4), 
and Al—)^A12 (arm-c5). The nomenclature here follows 


Watt & Mundy (1999) and Liu et al. (2012b), if the 
referred structures have been detected. The core Al is 
harboring OB stars, which create an UG Hll region. We 
refer to Liu et al. (2012b) for a preliminary study of the 
gas mass in the <0.1 pc scale structures. The molecu¬ 
lar arm-c4 appears spiraling out and connected with the 
parsec scale molecular arm-Sl and S2. The molecular 
arm-c2 may be connected with the previously observed 
molecular arm-N in the north (Liu et al. 2012b), which 
cannot be imaged well with our ALMA 12m Array data 
due to the smaller primary beam size. Following the con¬ 
vention in Liu et al. (2012b), we refer to the OB stars 
that are embedded in either core Al or core A2 as the 
central OB cluster hereafter. For the rest of localized 
intermediate- or high-mass (proto)stars and their par¬ 
ent molecular cores in G33.92-b0.11 A, we refer to them 
as the satellite high-mass stars and satellite cores. The 
previously identified core A4 is now resolved into two 
components A4n and A4s. In addition to the aforemen¬ 
tioned gas structures, we also resolved the ^0.1 pc scale 
gas mini-arms connecting with the cores. The gas mini¬ 
arms connecting A2 from its west (mini-arm-A2w) have 
been resolved in the previous SMA observations (Liu et 
al. 2012b). There are mini-arms connecting Al from its 
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Fig. 7. — The intensity weighted velocity dispersion maps of the selected molecular lines (color), overlaid on the 1.3 mm dust continuum 
image (contours). The velocity dispersion can be converted to full width of half maximum (FWHM) by multiplying the values by a factor 
of ~2.35. Contours are the same as those in Figure |4] We note that the blending of the CH 3 CN K=0, 1 components may artificially 
enhance the observed velocity dispersion in some regions. The position-velocity slices used for kinematic analysis are lab elled wi th bl ack 
and blue lines in the top middle and top right panels. The generated position-velocity diagrams will be presented in Figure fTOl 1131 and HU 
The colors of these labels are only for avoiding confusion. 


north and northwest (mini-arm-Ain, mini-arm-Alnw), 
which can only be marginally seen in Figure [T] owing to 
blending with other structures. We will address more 
about these gas mini-arms based on the spectral line re¬ 
sults (Section l3.2p . 

The high dynamic range of the ALMA observations 
confirms that the two highest mass cores Al and A 2 in 
the center are not round. The core Al appears slightly 
elongated in the southeast-northwest direction. This 
may imply that it is an inclined flattened structure, or 
there possesses eccentric gas orbits. But we cannot rule 
out the possibility that there is internal fragmentation in 
core Al, or that there are more than two cores blended 
in the projected line-of-sight. We cannot resolve the in¬ 
ternal structure of all detected cores. Some, or all of the 
individual cores may eventually form a cluster of stars. 
The complex morphology resolved by ALMA has sug¬ 
gested that the accretion flows in the massive molecular 
clump is highly dynamical/chaotic. In fact, the resolved 
arm-like features may resemble the scaled-down version 
of eccentric gas arms/arcs connected with the 2-4 pc scale 


Galactic circumnuclear disk/ring (Liu et al. 2012 d). 

3.2. Spectral lines 

We present the full spectral window spectra generated 
from the central 10" field of our ALMA 12 m-f-ACA ob¬ 
servations, in Figure [3l The identified line species are 
listed in Table [TJ In this work, we focus on the discus¬ 
sion of the CH 3 CN J=12-ll K-ladders, the 5-4 line, 
and the DCN 3-2 line. The velocity integrated intensity 
maps (i.e. moment 0) of these lines are shown in Fig¬ 
ure |H Figure [5] provides the side-by-side comparison of 
the blowed-up line images and the 1.3 mm dust contin¬ 
uum image in the central parsec scale area. The intensity 
weighted velocity maps (i.e. moment 1 ), and the inten¬ 
sity weighted velocity dispersion maps (i.e. moment 2 ) of 
these lines are presented in Figures [6H71 We present the 
velocity channel maps of the DCN 3-2 line, and of the 
^^CS 5-4 line in Figure [ 8 ] and [9l The position-velocity 
(PV) diagrams of these two lines, generated from the se¬ 
lected slices (see Figure [7]), are shown in Figure [iQl We 
present the velocity channel maps of the lowest and the 
highest transitions of the detected CH 3 CN J= 12 -ll K- 
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Fig. 8. — The velocity channel maps of the DCN 3-2 line (color), overlaid on the 1.3 mm dust continuum image (contour). Contours are 
1.5 mJy beam“^ X [1, 16]. 


ladders in Figures [TT] and [T2j The PV diagrams for these 
CH 3 CN transitions are provided in Figure [13] and Figure 
m We introduce the observed intensi ty dist rib ution a nd 
velocity profiles separately in Section 13.2.11 and 13.2.21 

3.2.1. Distribution 

The spatial distribution of the dense gas tracer DCN 
3-2 is in excellent agreement with the distribution of the 
1.3 mm dust emission structures (Figure IH [ 8 )). Strong 
DCN 3-2 emission is seen on the previously identified 
molecular clump C (Liu et al. 2012 b), all aforemen¬ 
tioned molecular arms, and shows several emission peaks 
coincide with the molecular cores A 1 -A 6 and All. The 
^^CS 5-4 line emission is also detected towards all these 
structures. However, the integrated ^^CS 5-4 map does 
not show enhanced emission from arm-c2 and arm-c3. 
From Figure m the most prominent ^^CS 5-4 emission 
structures are the mini-arm-Aln, mini-arm-Alnw, arm- 
c5, mini-arm-A2w, and the north part of arm-cl. The 
former three (mini-)arms engulf the converging point be¬ 
tween arm-c4 and core Al, and the later two surround the 
core A 2 . The ^^CS 5-4 enhanced mini-arms surrounding 
core Al are all redshifted relative to the systemic velocity 
107.6 kms“^ (more about the systemic velocity 
see Liu et al. 2012b and references therein). They start 
to emerge from the ambient gas structures in the veloc¬ 
ity channel of 108.28 kms“^ (Figure |9|), and becomes 
the brightest structures in the velocity channels redder 
than 110.32 kms“^. The redshifted mini-arm-Alnw can 


directly be seen in the average velocity map (Figure | 6 |), 
which artificially present a high velocity dispersion be¬ 
cause of blending with the ambient gas structures along 
the line-of-sight (Figure [7j). Toward the same area, we 
also found high NH 3 rotational temperature (Liu et al. 
2012 b). The mini-arm-Aln and mini-arm-Alnw are also 
weakly detected in the velocity maps of DCN 3-2 (Figure 
ini Figure [ 8 ]), which may indicate that these warm mini¬ 
arms are the result of the quickly dynamically processed 
gas streams. The redshifted motions of the mini-arm- 
Alnw can only be seen in its west end in the average 
velocity map of DCN 3-2 (Figure | 6 l top left), because of 
large confusion with the bulk motions of the ambient gas. 
The projected geometry of core Al together with arm-c4, 
mini-arm-Alw and mini-arm-Alnw highlighted by ^^CS 
5-4, resembles the scaled-down geometry of C33.92-f-0.ll 
A connected with the arm-N, SI, and S 2 seen in the pre¬ 
vious VLA NH 3 images and the SMA images of the CO 
2-1 isotopologues (Liu et al. 2012b). The ^^CS 5-4 line 
does not directly trace core A5, but some extended gas 
structures adjacent to it. This may suggest that core A5 
is surrounded by extended shocked gas. 

We only detected the K=0, 1, 2, 3 components of the 
CH 3 CN J= 12 -ll lines with upper level energies of 69 K, 
76 K, 97 K, and 133 K, respectively. These lines are en¬ 
hanced towards core Al, A 2 , and A5, but an extended 
and fainter component is seen in the entire ^ 0.6 pc area 
in C33.92-f0.ll A. In particular, all K components of 
CH 3 CN J= 12 -ll trace emission in mini-arm-A2w with 
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Fig. 9. — The velocity channel maps of the 5-4 line (color), overlaid on the 1.3 mm dust continuum image (contour). Contours are 

1.5 mJy beam“^ X [1, 16]. 


an extension to its northwest. The color composited im¬ 
ages in Figure [5] clearly demonstrate the similarity and 
the difference between the presented molecular gas trac¬ 
ers. 

3.2.2. Velocity Profile 

DCN trace dense gas, where the ices have been recently 
evaporated and injected into the gas phase. Since the 
destruction of deuterated neutral species in warm gas 
is slow (^10^-10^ yrs), DCN represents a fossil record 
of the deuterium content in the ices when the region 
was cold (Rodgers & Millar 1996). The 5-4 line 

traces large gas volume densities and large column den¬ 
sities, due to its high critical density and low abundance. 
Therefore, these two lines can probe the motions in the 
large scale gas ffows, without being seriously confused by 
localized (proto)stellar feedback. 

The velocity maps of the DCN 3-2 line and the ^^CS 
5-4 line consistently trace a southeast-northwest velocity 
gradient in the ^0.6 pc scale massive molecular clump 
G33.92-f0.llA (Figure [6l[Hip). The parsec scale molecu¬ 
lar arm-Sl,S2 appear blueshifted in the south and south¬ 
east ends, and become redshifted when moving closer 
to the dense molecular clump. Spatially, this is the most 
easily seen in the ^^CS 5-4 velocity channel map. The gas 
motion within the molecular arms are presented in the 
position-velocity (PV) diagrams of the ^^CS 5-4 and the 
DCN 3-2 lines (Figure HO]). We carefully placed the PV 
cuts, such that we can see the velocity gradient within the 


targeted molecular arms, without being confused much 
by the ambient dense gas structures. The locations of the 
presented PV cuts are drawn in the top middle panel of 
Figure 0 The PV diagram from the ‘cut c4d’ (Figure 0 
m is generated from the high spectral resolution (0.16 
kms“^) ^^CS 5-4 data taken with the ALMA 7m Array 
alone (6>majx^min=7".lx4".5, P.A.=89°), which is more 
sensitive to the mildly blueshifted motion in the south 
because of the larger field of view than the ALMA 12m 
Array observations. 

From Figure [TOl we see that although the DCN 3-2 
and the ^^CS 5-4 lines do not trace the same localized 
structures, the large scale velocity gradients traced by 
these two lines appear consistent. On the small scale, we 
see that the ^^CS 5-4 line traces slightly more redshifted 
emission by ^1 kms“^ in the molecular gas mini-arm- 
Alnw, as shown by the PV ‘cut Alw’. 

From ‘cut cF, ‘cut c2’, ‘cut c4a’, and ‘cut c5’, we do not 
resolve the Keplerian-like velocity profile. Rather than 
converging towards the cloud systemic velocity 
outer radii, we see the approximately linear velocity gra¬ 
dients in the molecular arms-cl, c2, and c5, which show 
the high velocity offsets at locations further away from 
the center of G33.92-f0.ll A. As a result, the average 
velocity map of 5-4 (Figure |6l top middle) trace 

the very blueshifted and redshifted motion in southeast 
and northwest of G33.92-f0.ll A, respectively. In the 
average velocity map of DGN 3-2 (Figure |6l top left), 
the redshifted motion in the northwest is confused by 
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Fig. 10. — Position-Velocity (PV) diagrams of the 5-4 (grayscale and magenta contours) and the DCN 3-2 lines (blue contours). 

Contours are 4.8 mjy beam“^(3cr)x [-1, 1, 2 , 4, 8, 16, 32, 64], except for the DCN 3-2 contours in the ‘cut Alw’ panel, which are 4.8 
mJy beam“^(3(7)x [1, 2, 3, 4, 5, 6]. Dashed lines label the n;s^=107.6 kms“^. The PV cuts are drawn in Figure [Jl with diamonds 
indicating the end points of the cuts (i.e. positive angular offset). Zero position offset is defined at the center of the presented PV cuts. 
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Fig. 11.— The velocity channel maps of the CH3CN J=12-ll, K=0 line (color), overlaid with the 1.3 m dust continuum image (contour). 
Contours are 1.5 mjy beam“^ x [1, 16]. 
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some diffused blueshifted high velocity features, which 
have been reported in the previous observations (Liu et 
al. 2012 b). From ‘cut c3a’ and ‘cut c3b’ (Figure [TO]), 
we see that the molecular arm-c3 is redshifted relatively 
to , and shows increasing velocity towards the north¬ 
west. The higher redshifted feature around the -0.8" off¬ 
set of ‘cut c3a’ is likely contributed by arm-c2. 

From the CH 3 CN velocity maps (Figure [ 6 ]), we see 
that the core Al is in general redshifted relative to core 
A 2 . Gaussian fits of the GH 3 GN spectral line profiles 
found the visr of cores Al and A 2 to be 108.93±1.19 


kms“^, and 107.55±0.83 kms“^, respectively. However, 
the GH 3 GN spectral profile in core Al is red-skewed, 
which may be due to spatial blending with the redshifted 
motion of mini-arm-Aln, or self-absorption. The cen¬ 
troid velocity of core Al determined by averaging the 
observed minimal and maximal terminal velocities, is 
^108.47 kms“^. The mini-arm-A2w has a blueshifted 
velocity of 107.19T0.83 kms“^. 

The comparison of the GH 3 GN and the DGN PV dia¬ 
grams show that the velocity offsets of core Al and A 2 
relative to are consistent with the global velocity gra- 
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dient within the 0.3 pc radius in G33.92+0.11 A (Figure 
m. These two massive cores are likely originated from 
the large scale accretion flow, and remain dynamically 
coupled to it. In addition, on the ^0.05 pc scale, we 
also observed the southeast-northwest velocity gradient 
within core Al. From the velocity channel maps of the 
K=0 component, we can see the emission peak in core 
Al moving from the southeast to the northwest, start¬ 
ing from the velocity channel of 106.18 kms“^ to the 
velocity channel of 108.43 kms“^ (Figure [12]). This mo¬ 
tion is blended with a more redshifted component, which 
is moving from south to north. The velocity channel 
map of the K=3 component trace better the south-north 
velocity gradient. We hypothesize that the south-north 
component is tracing the inner part of the mini-arm-Aln, 
although we cannot rule out the possibility that there is 
another redshifted core. The ‘cut Al’ panels in Figure 
[Ml show the PV diagrams of the K=0, 1, and 3 com¬ 
ponents of the CH 3 CN J= 12 -ll line, which were gen¬ 
erated from the PV cut closely following the southeast- 
northwest component in core Al (see also Figure 0 top 
right panel). We also do not resolve the Keplerian veloc¬ 
ity profile in core Al from the CH 3 CN PV diagrams (Fig¬ 
ure [M]). However, the highly accelerated gas in the center 
of core Al may already be photoionized, and therefore 
cannot be detected in CH 3 CN emission. The absence 
of high velocity CH 3 CN emission may also because of 
chemical segregation (see also Jimenez-Serra et al. 2012 ). 
Away from the center of core Al, CH 3 CN may not be ex¬ 
cited or is chemically less enhanced due to the lower gas 
temperature. 

We tentatively interpret the marginally resolved 
southeast-northwest velocity gradient by the rotational 
motion in Al. If this is the case, the consistently de¬ 
tected southeast-northwest velocity gradients within core 
Al and on the larger scale will be strongly indicative to 
a coherent gas accretion flow from large to small scale 
(more in Section lT3|) . These hypotheses can be tested 
in future higher angular resolution observations. We do 
not have enough sensitivity and angular resolution at the 
moment to resolve the velocity gradient within core A 2 
(and core A5). 


4. DISCUSSION 

We derive gas mass distribution based on the dust con¬ 
tinuum images, in Section 14.11 In Section 14.21 we dis¬ 
cuss the stability of the dense gas cores, which are likely 
the current or future sites of intermediate- or high-mass 
star-formation. Our interpretation of the resolved ve¬ 
locity field is provided in Section 14.41 which is based on 
our knowledge of gas mass distribution. We discuss the 
physical implication of our observational results in Sec¬ 
tion [441 


4.1. Matter Distribution 

The gas mass as well as column density in the ex¬ 
tended structures can be approximated by the optically 
thin dust emission formula 


Mh2 


2\^RapD‘^ 

3hcQ{X)J{X,Td) ^ 


( 1 ) 


where R is the gas-to-dust mass ratio, a is the mean grain 
radius, p is the mean grain density, D is the distance to 


the target, Q{X)(xX~^ is the grain emissivity, is the 
dust temperature, S{X) is the flux of the dust emission at 
the given wavelength, J{X,Td) = l/[ex.-p{hc/XkBTd) — 1 ] 
(Hildebrand 1983). The c, h, and ks are the light speed, 
the Planck constant, and the Boltzmann constant, re¬ 
spectively. Following Lis at al. (1998), we adopt the stan¬ 
dard values i?= 100 , a= 0.1 /im, p=3 gcm“^, (5(A=350 
pm)=l X10“"^. Based on the previous extensive NH 3 sur¬ 
vey toward the high-mass star-forming regions (Lu et 
al. 2014), and our previous NH 3 rotational temperature 
measurements on this region (Liu et al. 2012b), we adopt 
a constant Td=20 K on the large scale, and Tk^30 K in 
the inner 0.3 pc (^9") radius of G33.92+0.11 A. 

The gas mass of G33-north and G33-south in the poly¬ 
gon regions in Figure [T] are 9700 Mq and 5300 Mq, 
respectively. The gas column density in these two re¬ 
gions are in the range of A^hi+H 2 ^ 2 . 5-7.1x10^^ cm 
with a mean of ^4.4x10^^ cm“^. Given the assumed 
gas temperature and column density, and assuming that 
the width of the filaments are ^ 1-2 pc, the correspond¬ 
ing thermal Jeans length and Jeans mass are 0.20-0.47 pc 
and 7.2-17 Mq, respectively; the local free-fall collapsing 
timescale tff for the filaments is 0.42-1.0 Myrs. With the 
resolution of SHARG 2 /GSO (^0.34 pc) an d sen sitivity 
( 3 cr mass sensitivity ^ 10.5 Mq, see Section HHi we did 
not convincingly resolve the regularly spaced massive gas 
clumps in G33-north and G33-south (e.g. as seen in the 
more evolved OB cluster-forming region GlO.6-0.4, see 
Liu et al. 2012a). Nevertheless, there are a few local¬ 
ized over intensities at Jeans mass scales. The gas struc¬ 
tures western of G33.92+0.11 A are embedded within 
the localized gas clumps but are difficult to identify with 
the limited angular resolution of the SHARG 2 /GSO im¬ 
age. The enclosed gas mass in the central ^5 pc ra¬ 
dius is ^O.J^g gXlO^ Mq, where the error incorporates 
the uncertainties in dust temperature/opacity, as well as 
the foreground/background confusion. The global free- 
fall collapsing timescale for the region inside a 5 pc ra¬ 
dius is ^ 1.1 Myrs. We note that the detected gas mass 
within the central 5 pc radius in G33.92-l-0.il is com¬ 
parable to that in the Galactic mini-starburst region 
W49N (Galvan-Madrid et al. 2013). Therefore, we think 
G33.92-f0.ll has the potential of forming a very lumi¬ 
nous OB cluster in the future. We also note that the 
central ^0.3 pc radius massive clump G33.92+0.11 A is 
comparably massive (>3500 Mq) as the ^5 pc scale fil¬ 
aments G33-north/south, which indicates the high gas 
concentration in the cloud center. 

Finally, we summarize this section with two resolved 
features in common, from large to small scales. First, 
we detect the brightest sources at the center of the large 
scale 350 pm map, and of the high angular resolution 
ALMA 1.3 mm image (Figures HJ | 2 j). Within the 5 pc 
radius, the 350 pm image shows that ^3.8% of the gas 
mass is concentrated within 0.4% of the projected area 
(i.e. G33.92+0.11 A). Within source A inside a 0.3 pc 
radius, ALMA resolved that ^5.7% of the gas mass is 
concentrated to 0.7% of the projected area (i.e. core 
Al). Such a (self-)similarity in fraction of gas mass con¬ 
centrated in projected area may infer an anchored grav¬ 
itational collapse from large to small scales. Although 
these massive gas structures (i.e. G33.92+0.11 A and 
core Al) do not dominate the gas mass in the referred 
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Fig. 13.— Position-Velocity (PV) diagram for the CH3CN J=12- 
11, K=0,l,3 lines (magenta contour and grayscale) made within 
the central 0.3 pc radius in G33.92+0.11 A, overlaid with the PV 
diagram of the DON 3-2 line (cyan contour) generated from the 
same PV cuts. Dashed lines label the rJ^s-r^lOT.G kms“^. The 
PV diagrams of CH3CN lines are generated from the Briggs ro- 
bust=0 weighted image cubes, which has a synthesized beam of 
6>majX^min=0"-73x0".51 (P.A.=77°). The K=0 and 1 compo¬ 
nents of CH3CN J=12-ll are presented in the same panel, of which 
the velocity axis is alined with that of the K=0 component. The 
PV cuts are drawn in Figure [T] with diamonds indicating the end 
points of the cuts (i.e. positive angular offset). Magenta contours 
are 10.5 mJybeam“^(3cr)x[-l, 1, 2, 3, 4, 5, 6, 7, 8]. Cyan contours 
are 4.8 mJy beam”^ (3(7) x [-1, 1, 2, 4, 8, 16, 32]. 
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Fig. 14.— Position-Velocity (PV) diagram for the CH3CN J=12- 
11, K=0,l,3 lines (magenta contour and grayscale) made in core 
Al. Dashed lines label the kms”^. The PV cuts 

are drawn in Figure [Jl with diamonds indicating the end points 
of the cuts (i.e. positive angular offset). Magenta contours 
are 10.5 mjy beam”^ (3(7) x[-l, 1, 2, 3, 4, 5, 6, 7, 8]. These 
PV diagrams of CH3CN lines are generated from the Briggs ro- 
bust=0 weighted image cubes, which has a synthesized beam of 
0majXVin=O".73xO".51 (P.A.=77“). 


derive one-dimensional virial velocity dispersions 
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areas, locally they can be viewed as the dominant com¬ 
pact gravitational source. Both G33.92+0.11 A and core 
Al have companions with comparably lower masses. Sec¬ 
ond, we resolved the parsec scale molecular gas arms at 
where the ^5 pc scale gas streams converge to the mas¬ 
sive molecular clump. We also resolved the molecular gas 
mini-arms, which connect the <0.1 pc dense molecular 
core with the parsec scale molecular arms. The length 
of the observed molecular gas arms and mini-arms are 
^2-3 times the size of G33.92-b0.11 A and core Al, re¬ 
spectively. 

4.2. Gravitational stability of the satellite eores 

To examine whether the observed satellite cores in 
molecular arms can collapse to form high-mass stars, we 


for cores A3, A4n, A4s, A5-6, and A9-13 based on the 
formulation in Williams et al. (1994), and compare with 
the line-of-sight velocity dispersions {ay) traced by DGN 
3-2. In Equation [21 the parameter a is the geometric 
factor which equals to unity for a uniform density profile 
and 5/3 for an inverse square profile, M is the gas mass, 
G is the gravitational constant, and Reff is the effective 
radius of the dense clump. We suggest that the DGN 
3-2 line may be a particularly good dense gas tracer for 
this purpose, since DGN can survive 10^“^ years after the 
onset of (proto)stellar feedback (Rodgers & Millar 1996). 
To measure core masses and sizes, we fit 2-dimensional 
Gaussians to the 1.3 mm dust continuum images (Figure 
ED using the GAS A task imfit, and then estimate gas 
mass based on Equation 1. We assume the dust temper¬ 
ature K, and adopt /3=1. To take into account of 

the effects of ambient gas contamination and the missing 






















16 


Liu et al. 


flux, we perform Gaussian fittings to both the ALMA 
12 m+ACA Array 1.3 mm dust continuum image, and 
the dust continuum image generated from the 12 m Ar¬ 
ray data alone. For each of the cores, we extract the DCN 
line spectrum in the area defined by the 2 -dimensional 
Gaussian fitting, and then perform 1 -dimensional Gaus¬ 
sian fitting to the extracted spectrum to measure (jy . The 
results are presented in Figure [161 

We found that the derived Vyir is in general larger than 
(7y. This may imply that (1) the core masses are signifi¬ 
cantly overestimated by a factor of ^2.5 or more, or ( 2 ) 
the gas motions do not provide sufficient support against 
the local gravitational collapse, although we cannot rule 
out the case (3): all cores are rotationally flattened and 
the dominant gas motions are perpendicular to the line- 
of-sight. Case ( 1 ) may be true if large dust grains are 
present in the satellite cores, such that locally the value of 
/3 is close to 0 . Another option is that actual gas-to-dust 
mass ratio is a few times lower than the assumed value 
(100). Our data do not allow discriminating these cases. 
The values of /3 in satellite cores can be measured by 
performing dust continuum observations at other wave¬ 
lengths (e.g. 0.85mm or 3 mm). Case (3) needs to be 
tested by future higher angular resolution observations. 
We tentatively favor case ( 2 ), or a combination of cases 
(2) and (3) for some cores, since there were localized high 
velocity CO 2-1 outflow features detected in the previ¬ 
ous observations (Liu et al. 2012b). In any case, the 
evidence makes unlikely that the supersonic gas motions 
in the satellite cores are much faster than the actual Vyiy , 
such that these cores are not stable. 

4.3. Velocity Fields 

The observed gas motions in G33.92-h0.11 A, in par¬ 
ticular, in arms-cl, c2, c4 and c5, may be explained by 
rotational and free-fall motions, accelerated by the en¬ 
closed matter distribution. For example, we consider 
that the enclosed mass M(r) within the radius r can be 
approximated by M(r) = Mq • r^. Surrounding the ^0.3 
pc scale flattened rotating massive molecular clump, the 
value of n may be in between 1-3 (c.f. the solution for 
a collapsing slowly rotating isothermal cloud given by 
Terebey et al. 1984). The gas velocity can be estimated 
by Based on the observed in¬ 

creasing velocity offset at larger radii, the value of n may 
be approximately in between 2-3, which is expected in 
the region where the large scale gas streams are converg¬ 
ing to the flattened rotating structure. The PV ‘cut c4a’ 
does not show a clear velocity gradient in angular off¬ 
sets of -3"-0".5 (Figure Ho]), which is likely because the 
orientation of this PV cut is perpendicular to this ve¬ 
locity gradient of the rotating clump (c.f. Figure [ 6 ] and 
III). At larger angular offset, the ‘cut c4a’ detects the ex¬ 
terior redshifted infall (more below). We note that this 
analysis is preliminary because we do not know exactly 
how the inclination of the individual molecular arms vary 
with position. Defining the matter distribution as M(r) 
on the <0.3 pc scale, and deprojecting the observed line- 
of-sight velocity, are not possible because of this reason. 
Due to the complex morphology, the matter distribution 
may also be poorly approximated by M(r). The center 
of mass of G33.92+0.11 A is likely located in between 
cores Al and A 2 , but cannot be precisely determined 


at this moment. In addition, our observations do not 
allow us to separate the rotational motion from infall. 
Nevertheless, qualitatively we think our interpretation is 
physically plausible and has to be considered in other 
case studies. 

On smaller scales, the observed projected velocity sur¬ 
rounding core Al is ^ ±1.00 kms“^ relative to its sys¬ 
temic velocity n/s^=108.47 kms“^, at the 6 . 8 ± 0 . 8 xl 0 “^ 
pc projected radius. Assuming that the embedded 100 - 
300 Mq gas mass (Liu et al. 2012 b) in the core Al is 
dominant, and assuming that the surrounding gas in the 
< 0.01 pc radius follows the nearly circular orbits, the 
rotating speed ^^^(0.0068 pc) can be estimated by vr = 
[(G ■ (100-300 M 0 ))/(O.OO 68 pc)]i/2 =7.7-13 kms-i. We 
therefore can interpret the southeast-northwest velocity 
gradient around Al by the rotation inclined by 4°-7° 
along the position angle of ^45°. The CH 3 CN velocity 
field surrounding the core A 2 is more complicated, in par¬ 
ticular, in the mini-arm-Aw. Nevertheless, we found that 
at the 0.035-0.05 pc position offset from core A 2 , the ob¬ 
served ~0.4 kms“^ velocity offset of the mini-arm-A 2 w 
from the centroid velocity of core A 2 , can be consistently 
interpreted by the gravitationally accelerated motion by 
core A 2 (and Al), projected to a similar inclination angle 
to that of core Al. 

If the velocity gradient on the >0.3 pc scale is consis¬ 
tent with that on the <0.3 pc as well as the southeast- 
northwest velocity gradient within core Al, we could ex¬ 
pect the entire molecular arm-Sl and S 2 (i.e. southern 
than deck (J2000)~00°55'23") to be blue shifted. We 
present the PV ‘cut c4b’, ‘cut c4c’, ‘cut c4d’ in Figure 
HQl which trace the gas motions in the molecular arm- 
Sl and S 2 . We instead found that the line-of-sight mo¬ 
tion of these two arms converge smoothly from the red- 
shifted velocity in the north to the blueshifted velocity 
in the south and southeast. These results are consistent 
with the previous 0.6 kms“^ velocity resolution observa¬ 
tions of the NH 3 (1,1) hyperflne inversion lines (Liu et al. 
2012b). In the marginally spatially resolved system, this 
approximately linear velocity gradient is most commonly 
interpreted by the rotational motion of a toroid. In the 
particular case of G33.92±0.11, we are motivated by the 
resolved matter distribution, and the morphology of the 
arms, to alternatively consider that the parsec scale ve¬ 
locity gradient in the north-south direction can be the 
free-fall dominant motion following the eccentric orbits. 
The resolved velocity gradient from the molecular arm-c4 
to arms-Sl, S 2 , is qualitatively similar to the numerical 
simulation for the infalling filament connected with the 
low-mass cluster-forming region in Bonnell et al.(2008). 
The southern molecular arm may be gravitationally ac¬ 
celerated by the massive core Al, and therefore turns 
into the redshifted velocity in regions close to Al. From 
both the average velocity map and the PV diagrams of 
5-4 (Figure (61 fTO]) . we observe the smaller velocity 
scale outside the 0.3 pc radius of G33.92±0.11 A, than 
the velocity of the gas arms interior to it. This can be 
understood if the observed gas structures within the 0.3 
pc radius of G33.92±0.11 A dominate the enclosed mass 
within the ^1 pc radius. 

We summarize our working hypothesis to interpret 
the resolved structures and kinematics in the schematics 
model in Figure HSl We think the molecular gas streams 
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are converging from the larger than 1 pc radius towards 
the center of G33.92+0.11 A, due to gravitational infall 
acceleration. The gas motions become dominated by the 
rotational motion on the ^0.3 pc scale. Within the 0.3 
pc scale radius, we detect the gas mini-arms connecting 
to the central core A 1 and A 2 , which are rotating about 
each other. These gas mini-arms show <1 kms“^ ve¬ 
locity offsets from the ambient dense gas, which can be 
explained by the gravitational acceleration and the tidal 
interaction of the dense core. However, we cannot yet 
rule out that the motion of the gas mini-arms is pow¬ 
ered by the expansional motion of the ionized gas, or 
the molecular wind/jet, which needs to be examined in 
higher resolution observations. 

4.4. Physical Implications 

The magnetic field in rapidly collapsing molecular 
clumps is presumably weak. The initial turbulence is 
rapidly dissipated in such dense regions, although can 
be replenished by the protostellar turbulence in a later 
evolutionary stage (Wang et al. 2010). The rotational 
motion can play an important role in supporting collaps¬ 
ing molecular clumps, which permits localized fragmen¬ 
tation and star-formation in longer timescales. It is not 
yet clear how to transport angular momentum outward 
efficiently on spatial scale of a fraction of a parsec. The 
excess of angular momentum may leave its footprint on 
the morphology of the molecular clumps, and may there¬ 
fore bias the upper end of the core/stellar mass function. 

Observationally, we have found clumpy rotating 
toroids and spiral-like interacting gas features toward 
the low-mass cluster-forming region LI287 (Juarez et 
al. in prep.), the OB cluster-forming molecular clumps 
(e.g. GlO.6-0.4: Liu et al. 2010; G35: Sanchez-Monge 
et al. 2013, Qiu et al. 2013; NGG6334V: Juarez et 
al. in prep.), and the Galactic mini-starburst region 
W49N (Welch et al.l987, Galvan-Madrid et al. 2013, 
etc). These results imply that this physical problem is 
spatially scalable, and the previously spatially resolved 
systems are the large-scale examples seen in a closer to 
face-on projection. However, toroids may also be the 
projected central OB cluster-forming cores surrounded 
by satellite cores. The critical spatial scales to look at, 
are the centrifugal radius (or radii), where the centrifugal 
force marginally balances the gravitational force from the 
enclosed molecular and stellar mass. We expect more of 
these systems to be discovered or be spatially resolved in 
high angular resolution ALMA observations. We hypoth¬ 
esize that the inward migration of satellite cores can ex¬ 
plain the formation of the extremely concentrated dense 
condensations, similar to what was resolved in NH 3 (3,3) 
satellite hyperfine line absorption against the UG Hll re¬ 
gion GlO.6-0.4 (Sollins & Ho 2005; Liu et al. 2010). 

Finally, we hypothesize that the spiraling arm-like 
asymmetry is essential in the gravitationally dominated 
accretion flow, and needs to be considered in modeling 
frameworks. These structures may either be induced at 
regions where the eccentric accretion streams collide with 
each other, or be created by tidal interactions. We note 
that tidal accretion of eccentric accretion streams, or 
tidally interacting hot cores, may deposit the hot molec¬ 
ular gas tracers from the hot cores to the entire parsec 
scale massive molecular clump. 

We would like to emphasize that the location of the 



Fig. 15. — A schematic model of the resolved region. Gray el¬ 
lipses are the projected iso-radius contours of an inclined circular 
rotating plane. The near and far sides of this inclined system are 
separated by the dashed line. The blue and red arrows show the 
spiraling ~0.5 pc scale molecular arms, and the mini-arms con¬ 
nected with cores Al and A2, respectively. The filled circles are 
the cores Al and A2, with color scales indicating the line-of-sight 
velocity gradients caused by their spinning motion. These two mas¬ 
sive cores are rotating about each other, such that the core Al is 
redshifted with respect to core A2. Observationally, the velocity 
distribution in A2 is confused with the blueshifted tidally interact¬ 
ing feature connecting from A2 to Al, which is indicated as the 
filled cyan area. We note that some of the arms can be off plane. 

objects in the region discussed favor the idea of the con¬ 
vergence of parsec scale gas accretion flows onto a central 
clump. Our high spectral resolution DGN 3-2 and ^^GS 
5-4 line observations suggest that this massive molecu¬ 
lar clump is fed by free-falling exterior gas streams (or 
arms). Based on the resolved images, we found that 
accretion flows in this region are highly dynamical and 
spatially non-uniform. Individual accretion gas streams 
likely carry different amount of specific angular momen¬ 
tum. Even in the case that the averaged specific angu¬ 
lar momentum on the large scale is negligibly small, it 
remains uncertain how the specific angular momentum 
carried by individual gas streams will alter the gas accre¬ 
tion and fragmentation. High angular resolution molec¬ 
ular line observations to study a large ensemble may be 
required to elucidate this issue. 
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Fig. 16.— Estimated one-dimensional virial velocity dispersions 
{vyir) and observed line-of-sight velocity dispersions (a y) o f satel¬ 
lite cores A3, 4n, 4s, 5, 6, 9, 10, 11, 12 and 13 fSection 14.2|) . Blue 
symbols are the measurements based on the ALMA 12m+ACA 
Arrays 1.3 mm dust continuum image, and gray symbols are the 
measurements based on the ALMA 12m Array 1.3 mm dust contin¬ 
uum image. Red dashed line shows Vyir=(Jv The horizontal error 
bars are ±0.34 kms“^ (i.e. half of the velocity channel spacing). 
The upper and lower bounds of the vertical error bars are derived 
by assuming the inverse-square radial density profile and the uni¬ 
form density distribution, respectively. We note that core 12 is 
seriously confused with the negative sidelobes in the 12m Array 
alone image, and therefore its flux cannot be measured. 


5. SUMMARY 

We observed the approximately face-on young 
OB cluster-forming region G33.92±0.11, using 
SHARC 2 /CSO, Herschel, and ALMA. On the >1 
pc scale, the SHARC 2 0.35 mm dust continuum image 
trace the filamentary gas streams converging to the cen¬ 
tral ^ 10 ^ Mq massive molecular gas clump G33.92±0.11 
A. The high resolution ALMA 1.3 mm dust continuum 
images reveal abundant dense satellite cores, which 
are closely associated with the ^0.5 pc scale molecular 
gas arms orbiting the two central highest mass cores. 
We found that the DGN 3-2 line emission correlates 
well with the 1.3 mm dust continuum emission, while 
the 5-4 line shows further enhanced emission 

towards the gas mini-arms connecting with the central 
massive cores and the potentially shock heated gas. 
These two lines consistently trace a southeast-northwest 
global velocity gradient within the 0.3 pc radius in 
G33.92±0.11 A, which may be interpreted by the 


rotational and infall motion. The small line width and 
velocity differences between dense features strongly 
supports the interpretation of a flattened object viewed 
nearly face-on. The spiral arms would indicate that 
the flattening is mostly likely due to rotation, although 
the full rotation speed is hard to resolve because of 
the face-on geometry. The GH 3 GN lines appear to 
trace the same sense of rotation to the <0.05 pc scale 
hot core regions, which indicated that there may be 
a coherent accretion flow continuing from large to 
small scale. The parsec scale molecular arm-Sl, S 2 
traced by the 1.3 mm dust continuum emission, and the 
molecular line emission, may be the infall gas streams 
feeding massive molecular clump G33.92±0.11 A. We 
highlight the resolved arm-like features connecting with 
the massive localized structures, and hypothesize that 
the arm-like morphology is essential in star and cluster 
forming regions of all masses. We derived higher virial 
velocity dispersions than observed velocity dispersions 
for the satellite molecular cores embedded in molecular 
arms, which suggest that these cores are gravitationally 
unstable, and may collapse to form high-mass stars. 


This paper makes use of the following ALMA data: 
ADS/JAO.ALMA 2012.1.00387.S. ALMA is a partner¬ 
ship of ESO (representing its member states), NSF 
(USA) and NINS (Japan), together with NRG (Ganada) 
and NSG and ASIAA (Taiwan), in cooperation with the 
Republic of Ghile. The Joint ALMA Observatory is op¬ 
erated by ESO, AUI/NRAO and NAOJ. This material 
is based upon work at the Galtech Submillimeter Ob¬ 
servatory, which is operated by the Galifornia Institute 
of Technology. In the present paper we discussed ob¬ 
servations performed with the ESA Herschel Space Ob¬ 
servatory (Pilbratt et al. 2010), in particular employing 
Herschel’s large telescope and powerful science payload 
to do photometry using the PACS (Poglitsch et al. 2010 ) 
and SPIRE (Griffin et al. 2010) instruments. Herschel 
is an ESA space observatory with science instruments 
provided by European-led Principal Investigator consor¬ 
tia and with important participation from NASA. This 
research made use of Montage, funded by the National 
Aeronautics and Space Administration’s Earth Science 
Technology Office, Gomputation Technologies Project, 
under Cooperative Agreement Number NCC5-626 be¬ 
tween NASA and the California Institute of Technology. 
Montage is maintained by the NASA/IPAC Infrared Sci¬ 
ence Archive. HBL thanks Simon Radford and other 
CSO staffs for assisting the operation of SHARC 2 . HBL 
thanks ASIAA for supporting series of related researches. 
HBL acknowledge funding under MOST grant number 
103-2119-M-001-010-MY2. I.J.-S. acknowledges funding 
from the People Programme (Marie Curie Actions) of 
the European Union’s Seventh Framework Programme 
(FP7/2007-2013) under REA grant agreement number 
PHF-GA-2011-301538. CRZ acknowledges support from 
program CONACYT CB2010-152160, Mexico. 
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APPENDIX 

A. IDENTIFIED LINES FROM THE ALMA 12M + ACA OBSERVATIONS 

We present the observed molecular lines in the central 10" region in G33.92+0.11 A, in Table [TJ 
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TABLE 1 

Molecular Lines observed within our Cycle- 1 ALMA data 


Species 

Transition 

Frequency 

(MHz) 

Eu 

(K) 

C-C3H2 

33,0—^22,1 

216278.76 

19.5 

CCD 

N=3-2, J=7/2-5/2, F=9/2-7/2 

216372.83 

20.8 


N=3-2, J=7/2-5/2, F=7/2-5/2^ 

216373.32 

20.8 


N=3-2, J=5/2-3/2, F=7/2-5/2‘> 

216428.32 

20.8 


N=3-2, J=5/2-3/2, F=3/2-l/2 

216428.76 

20.8 

CH3CHO 

111,10 ^10i,9 E 

216581.94 

64.9 


111,10 ^10i,9 A 

216630.22 

64.8 

HDCS 

7o,7 ^ 60,6 

216662.43 

41.6 


72,6 ^62,5 

216931.37 

77.6 


^ 2,5 ^62,4 

217263.69 

77.6 

H2S 

22,0 ^2iq 

216710.44 

84.0 

CH3OH 

5i,4 ^42,2 

216945.60 

55.9 

SiO 

5^4 

217104.98 

31.3 

DCN 

3^2 

217238.63 

20.9 

13CN 

N=2-l, J=3/2-l/2, Fl=l-1, F=2-l 

217072.80 

15.7 


N=2-l, J=3/2-l/2, Fl=l-1, F=2-2 

217074.24 

15.7 


N=2-l, J=3/2-l/2, Fl=l-0, F=0-1 

217264.64 

15.7 


N=2-l, J=3/2-l/2, Fl=l-0, F=l-1 

217277.68 

15.7 


N=2-l, J=5/2-3/2, Fl=2-1, F=2-2 

217286.80 

15.7 


N=2-l, J=5/2-3/2, Fl=2-1, F=l-1 

217290.82 

15.7 


N=2-l, J=5/2-3/2, Fl=2-1, F=l-0 

217296.60 

15.7 


N=2-l, J=3/2-3/2, Fl=2-2, F=2-2 

217298.94 

15.7 


N=2-l, J=5/2-3/2, Fl=2-1, F=2-l 

217301.18 

15.7 


N=2-l, J=5/2-3/2, Fl=2-1, F=3-2 

217303.19 

15.7 


N=2-l, J=3/2-l/2, Fl=l-0, F=2-l 

217304.93 

15.7 


N=2-l, J=3/2-3/2, Fl=2-2, F=3-3 

217306.12 

15.7 


N=2-l, J=5/2-3/2, Fl=2-1, F=3-2 

217428.56 

15.7 


N=2-l, J=5/2-3/2, Fl=2-1, F=2-l 

217436.35 

15.7 


N=2-l, J=5/2-3/2, Fl=2-1, F=2-2 

217437.70 

15.7 


N=2-l, J=3/2-l/2, Fl=2-1, F=2-2 

217437.82 

15.7 


N=2-l, J=3/2-l/2, Fl=2-1, F=l-1 

217443.72 

15.7 


N=2-l, J=5/2-3/2, Fl=3-2, F=4-3 

217467.15 

15.7 


N=2-l, J=5/2-3/2, Fl=3-2, F=2-l 

217469.15 

15.7 


N=2-l, J=5/2-3/2, Fl=3-2, F=2-2 

217480.56 

15.7 


N=2-l, J=5/2-3/2, Fl=3-2, F=3-3 

217483.61 

15.7 

t-C2H50H 

53,3 ^42,2 

217803.69 

23.9 

c-HCCCH 

61,6 ^5 o, 5^ 

217822.15 

38.6 


5i,4 ^42,3 

217940.05 

35.4 

CH3CN 

12^11 K=4 

220679.29 

183.1 


12^11 K=3 

220709.02 

133.2 


12^11 K=2 

220730.26 

97.4 


12^11 K=1 

220743.01 

76.0 


12^11 K=0 

220747.26 

68.9 

CCS 

19^18 

231060.99 

110.9 

i^CS 

5^4 

231220.69 

33.3 

Ha 

n=30 

231900.93 

0.0 

Hea 

n=30 

231995.43 

0.0 

CH3OCH3 

13o,13 ^12i,i2 EE*^ 

231987.82 

80.9 

H2C34s 

7i,7 ^ 61,6 

232754.71 

57.9 


TABLE 2 
* 

“Blended with transition CCD N=3-2, J=7/2-5/2, E=5/2-3/2. 

^Blended with transition CCD N=3-2, J=5/2-3/2, F=5/2-3/2. 

“Blended with transition c-HCCCH 6o,6 

^The CH3OCH3 13o, 13 transitions from the AA, AE and EA species are blended 

with this line (frequencies are within ±0.2 MHz). 






